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ABSTRACT 

The  Hubble  Ultra  Deep  Field  (HUDF)  contains  a  significant  number  of  B-,  V-,  and  /'-band  dropout  objects, 
many  of  which  were  recently  confirmed  to  be  young  star-forming  galaxies  at  z  —  4-6.  These  galaxies  are  too 
faint  individually  to  accurately  measure  their  radial  surface-brightness  profiles.  Their  average  light  profiles  are 
potentially  of  great  interest,  since  they  may  contain  clues  to  the  time  since  the  onset  of  significant  galaxy  assembly. 

We  separately  co-add  V-,  /'-,  and  z'-band  HUDF  images  of  sets  of  z  —  4,  5,  and  6  objects,  pre-selected  to  have 
nearly  identical  compact  sizes  and  the  roundest  shapes.  From  these  stacked  images,  we  are  able  to  study  the 
average(d)  radial  structure  of  these  objects  at  much  higher  signal-to-noise  ratio  than  possible  for  an  individual 
faint  object.  Here,  we  explore  the  reliability  and  usefulness  of  a  stacking  technique  of  compact  objects  at  z  —  4-6 
in  the  HUDF.  Our  results  are:  (1)  image  stacking  provides  reliable  and  reproducible  average  surface-brightness 
profiles;  (2)  the  shape  of  the  average  surface-brightness  profile  shows  that  even  the  faintest  z  —  4-6  objects 
are  resolved ;  and  (3)  if  late-type  galaxies  dominate  the  population  of  galaxies  at  z  —  4-6,  as  previous  Hubble 
Space  Telescope  studies  have  shown  for  z  <  4,  then  limits  to  dynamical  age  estimates  for  these  galaxies  from 
their  profile  shapes  are  comparable  with  the  spectral  energy  distribution  ages  obtained  from  the  broadband  colors. 

We  also  present  accurate  measurements  of  the  sky  background  in  the  HUDF  and  its  associated  I  n  uncertainties. 

Key  words:  galaxies:  formation  -  galaxies:  high-redshift  -  galaxies:  structure 


1.  INTRODUCTION 

In  the  last  decade,  ground-  and  space-based  observations 
of  high-redshift  galaxies  have  begun  to  outline  the  process  of 
galaxy  assembly.  The  details  of  that  process  at  high  redshifts, 
however,  remain  poorly  constrained.  There  is  increasing  support 
for  the  model  of  galaxy  formation,  in  which  the  most  massive 
galaxies  assemble  earlier  than  their  less  massive  counterparts 
(e.g.  Cowie  et  al.  1996;  Guzman  et  al.  1997;  Kodama  et  al.  2004; 
McCarthy  2004).  A  detailed  analysis  of  the  “fossil  record”  of  the 
current  stellar  populations  in  nearby  galaxies  selected  from  the 
Sloan  Digital  Sky  Survey  (SDSS;  York  et  al.  2000)  provides 
strong  evidence  for  this  downsizing  picture  (Heavens  et  al. 
2004;  Panter  et  al.  2007).  The  increasing  number  of  luminous 
galaxies  spectroscopically  confirmed  to  be  at  z  —  6.5  (e.g. 
Hu  et  al.  2002;  Kodaira  et  al.  2003;  Kurk  et  al.  2004;  Rhoads 
et  al.  2004;  Stern  et  al.  2005;  Taniguchi  et  al.  2005),  or 
<0.9  Gyr  after  the  Big  Bang,  also  supports  this  general  picture. 
In  an  alternative  hierarchical  scenario,  arguments  have  been 
made  that  a  significant  number  of  low  luminosity  dwarf  galaxies 
were  present  at  these  times,  and  were  the  main  contributor  to 
finish  the  process  of  reionization  of  the  intergalactic  medium 
(Yan  &  Windhorst  2004a,  2004b).  However,  there  is  currently 
little  information  on  the  dynamical  structure  of  these  or  other 
galaxies  at  z  —  6.  It  is  not  clear  whether  these  objects  represent 
isolated  disk  systems,  or  collapsing  spheroids,  mergers  or  other 
dynamically  young  objects. 

Ravindranath  et  al.  (2006)  used  deep,  multi-wavelength 
images  obtained  with  the  Hubble  Space  Telescope  ( HST ) 
Advanced  Camera  for  Surveys  (ACS)  as  part  of  the  Great 
Observatories  Origins  Deep  Survey  (GOODS)  to  analyze  2D 
surface-brightness  distributions  of  the  brightest  Lyman  Break 
Galaxies  (LBGs)  at  2.5  <  z  <  5.  They  distinguish  various 


morphologies  based  on  the  Sersic  index  n,  which  measures 
the  shape  of  the  azimuthally  averaged  surface-brightness  pro¬ 
file  (where  n  =  1  for  exponential  disks  and  n  =  4  for  a  de 
Vaucouleurs  law).  Ravindranath  et  al.  (2006)  find  that  40%  of 
the  LBGs  have  light  profiles  close  to  exponential,  as  seen  for 
disk  galaxies,  and  only  ~30%  have  high  n,  as  seen  in  nearby 
spheroids.  They  also  find  a  significant  fraction  (~30%)  of  galax¬ 
ies  with  light  profiles  shallower  than  exponential,  which  ap¬ 
pear  to  have  multiple  cores  or  disturbed  morphologies,  sug¬ 
gestive  of  close  pairs  or  on-going  galaxy  mergers.  Distinction 
between  these  possible  morphologies  and,  therefore,  a  better 
estimate  of  the  formation  redshifts  of  the  systems  observed 
at  z  —  4-6  in  particular,  is  important  for  testing  the  galaxy 
assembly  picture,  and  for  the  refinement  of  galaxy  formation 
models. 

One  possible  technique  involves  the  radial  surface-brightness 
profiles  of  the  most  massive  objects — those  that  will  likely 
evolve  to  become  the  massive  elliptical  galaxies,  which  we  see 
in  place  at  redshifts  z  <  2  (Driver  et  al.  1998;  van  Dokkum 
et  al.  2003,  2004).  This  can  be  analytically  understood  in  the 
context  of  the  Lynden-Bell  (1967)  relaxation  formalism  and  the 
numerical  galaxy  formation  simulations  of  van  Albada  (1982), 
which  describe  collisionless  collapse  and  violent  relaxation  as 
the  formation  mechanism  for  elliptical  galaxies.  As  the  time 
scale  for  relaxation  is  shorter  in  the  inner  than  in  the  outer  parts 
of  a  galaxy,  convergence  toward  an  r<"  profile  will  proceed 
from  the  inside  to  progressively  larger  radii  at  later  times. 
Moreover,  Kormendy  (1977)  has  shown  that  tidal  perturbations 
due  to  neighbors  can  cause  the  radial  surface-brightness  profile 
to  deviate  from  a  pure  de  Vaucouleurs  profile  in  the  outer  parts  of 
a  galaxy.  This  implies  that  the  radius  where  surface-brightness 
profiles  start  to  deviate  significantly  from  an  rl/'n  profile  might 
serve  as  a  “virial  clock”  that  traces  the  time  since  the  onset  of 
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the  last  major  merger,  accretion  events,  or  global  starburst  in 
these  objects. 

Image  stacking  methods  have  been  used  extensively  on  X-ray 
(Brandt  et  al.  2001;  Nandra  et  al.  2002)  and  radio  (Georgakakis 
et  al.  2003;  White  et  al.  2007)  data  to  study  the  mean  properties 
(e.g.  flux,  luminosity)  of  well-defined  samples  of  sources  that 
are  otherwise  too  faint  to  be  detected  individually.  Pascarelle 
et  al.  (1996)  applied  such  a  stacking  method  to  a  large  number 
of  optically  very  faint,  compact  objects  at  z  =  2.39  to  trace  their 
“average”  structure.  This  approach  was  also  applied  by  Zibetti 
et  al.  (2004)  to  detect  the  presence  of  faint  stellar  halos  around 
disk  galaxies  selected  from  the  SDSS.  An  attempt  to  apply  this 
technique  to  high  redshift  galaxies  in  the  Hubble  Deep  Field 
(HDF;  Williams  et  al.  1996)  was  not  conclusive  (H.  Ferguson, 
private  communication)  due  to  the  poorer  spatial  sampling  and 
shallower  depth  of  the  HDF  compared  to  the  Hubble  Ultra  Deep 
Field  (HUDF;  Beckwith  et  al.  2006). 

In  this  paper,  we  use  the  exceptional  depth  and  fine  spatial 
sampling  of  the  HUDF  to  study  the  potential  of  this  image 
stacking  technique,  and  will  estimate  limits  to  dynamical  ages 
of  faint,  young  galaxies  at  z  —  4-6.  The  HUDF  reaches  ~1.5 
mag  deeper  than  the  equivalent  HDF  exposure  in  the  ('-band,  and 
has  better  spatial  sampling  than  the  HDF.  The  HUDF  depth  also 
allows  us  to  characterize  the  sky  background  very  accurately, 
which  is  critical  for  successfully  using  a  stacking  method  to 
measure  the  mean  surface-brightness  profiles  for  these  faint 
young  galaxies. 

This  paper  is  organized  as  follows.  In  Section  2  we  summarize 
the  HUDF  observations,  and  in  Section  3  we  discuss  the 
selection  of  our  z  —  4,  5,  and  6  samples.  In  Section  4  we 
describe  our  data  analysis,  which  includes  accurately  measuring 
the  1  a  sky-subtraction  error,  the  image  stacking  method  to 
generate  mean  surface-brightness  profiles,  and  our  test  of  its 
reliability.  In  Section  5  we  present  and  discuss  our  results  in 
terms  of  the  average  surface-brightness  profiles  of  z  ~  4-6 
galaxies,  and  in  Section  6  we  conclude  with  a  summary  of  our 
results. 

Throughout  this  paper  we  refer  to  the  HST/ACS  F435W, 
F606W,  F775W,  and  F850LP  filters  as  the  B-,  V-,  i'-,  and 
z'-bands,  respectively.  We  assume  a  Wilkinson  Microwave 
Anisotropy  Probe  ( WMAP )  cosmology  of  Qm  =  0.24,  Q.A  = 
0.76,  and  Hq  =  73  km  s-1  Mpc-1 ,  in  accord  with  the  most  recent 
3-year  WMAP  results  of  Spergel  et  al.  (2007).  This  implies  a 
current  age  for  the  Universe  of  13.65  Gyr.  All  magnitudes  are 
given  in  the  AB  system  (Oke  &  Gunn  1983). 

2.  OBSERVATIONS 

The  HUDF  contains  >100  objects  that  are  ('-band  dropouts, 
making  them  candidates  for  galaxies  at  z  —  6  (Bouwens  et  al. 
2004,  2006;  Bunker  et  al.  2004;  Yan  &  Windhorst  2004b). 
Similarly,  there  are  larger  numbers  of  objects  in  the  HUDF 
that  are  B-band  dropouts  (415  in  total),  or  V-band  dropouts 
(265  in  total),  and  are  candidates  for  galaxies  at  z  —  4  and 
z  —  5,  respectively.  Beckwith  et  al.  (2006)  and  Bouwens  et 
al.  (2007)  find  similar  number  of  B-  and  V-band  dropouts 
in  the  HUDF.  A  significant  fraction  of  these  objects  to  AB 
<  27  mag  has  recently  been  spectroscopically  confirmed  to 
have  redshifts  z  ~  4-6  through  the  detection  of  Lya  emission, 
or  identifying  their  Lyman  break  (Malhotra  et  al.  2005;  Dow- 
Hygelund  et  al.  2007).  We  discuss  our  detailed  dropout  selection 
criteria  below.  Despite  the  depth  (AB  <  29.5  mag)  of  the 
HUDF  images,  however,  these  objects  appear  very  faint  and 
have  little,  if  any,  discernible  structural  detail.  Visual  inspection 


Figure  1.  The  HUDF  number  counts  for  all  z  —  4,  5,  and  6  objects  before 
the  sub-selection  of  compact  isolated  z  —  4,  5,  and  6  objects  was  made.  The 
vertical  dotted  line  shows  the  magnitude  to  which  the  number  counts  of  all  these 
redshifts  are  complete.  The  area  of  the  HUDF  is  3.15  x  10-3  deg2. 

of  all  these  objects  shows  their  morphologies  to  divide  into 
four  broad  categories:  symmetric,  compact,  elongated,  and 
amorphous. 

3.  SAMPLE  SELECTION 

We  construct  three  separate  catalogs  for  these  z  —  4,  5,  and 
6  galaxy  candidates,  selecting  only  the  isolated,  compact,  and 
symmetric  galaxies.  We  exclude  objects  with  obvious  nearby 
neighbors,  to  avoid  a  bias  due  to  dynamically  disturbed  objects 
and  complications  due  to  chance  superpositions.  Figure  1 
demonstrates  that  our  completeness  limit  for  z  —  4  and  z  ~  5 
objects  is  AB  <  29.3  mag,  and  for  z  —  6  objects  it  is  AB  < 
29.0  mag.  Therefore,  all  three  catalogs  are  complete  to  AB  < 
29.0  mag,  which  is  equivalent  to  at  least  a  10cr  detection  for 
objects  that  are  nearly  point  sources.  For  each  object  in  our 
z  —  4,  5,  and  6  samples,  we  extracted  51x51  pixel  postage 
stamps  (which  at  (X'03  pix-1  span  1"53  on  a  side)  from  the 
HUDF  V-,  ('-,  and  z'-band  images,  respectively.  Each  postage 
stamp  was  extracted  from  the  full  HUDF,  such  that  the  centroid 
of  an  object  (usually  coincident  with  the  brightest  pixel)  was  at 
the  center  of  that  stamp. 

3.1.  The  z  —  4  and  z  —  5  Objects  (B-,  V  -Band  Dropouts) 

We  used  the  ('-band  selected  BVi'z!  HUDF  catalog 
(Beckwith  et  al.  2006)  to  select  the  z  ~  4  and  z  —  5  objects. 
With  the  HyperZ  code  (Bolzonella  et  al.  2000),  we  computed 
photometric  redshift  estimates,  using  the  magnitudes  and  asso¬ 
ciated  uncertainties  tabulated  in  the  HUDF  catalog.  All  objects 
with  3.5  s;  Zphot  ^  4.5  were  assigned  to  the  bin  of  z  —  4  candi¬ 
dates,  and  all  objects  with  4.5  ^  zphot  ^  5.5  to  the  bin  of  z  ~  5 
candidates. 

We  then  applied  color  criteria,  similar  to  those  adopted  by 
Giavalisco  et  al.  (2004),  to  select  the  B  (z  —  4)  and  V  (z  ~  5) 
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dropout  samples.  For  B -band  dropouts,  we  require 
(B  -  V)  Js  1.2  +  1.4  x  (V  -  z')  mag 
and  (B  —  V)  >  1.2  mag 
and  ( V  -  z')  ^1.2  mag. 

For  F-band  dropouts,  the  following  color  selection  was  applied: 

I(V  —  /')  >1.5  +0.9  x  (/'  —  z')  or  ( V  —  /')  >  2.0  mag 
and  ( V  —  /')  >  1.2  mag 
and  ( /'  —  z')  ^1.3  mag. 

We  note  that  only  objects  satisfying  both  color  and  photometric 
redshift  criteria  were  selected  in  our  samples.  Vanzella  et  al. 
(2006)  using  VFT/FORS2  observed  ~  1 00  B-,  V-.  and  /'-band 
dropout  objects  in  the  Chandra  Deep  Field  South  (CDFS) 
selected  based  on  above-mentioned  color  criteria  (Giavalisco 
et  al.  2004).  They  have  spectroscopically  confirmed  >90% 
of  their  high-redshift  galaxy  candidates.  Therefore,  we  expect 
only  a  small  number  (<10%)  of  contaminants  in  our  sample  of 
dropouts.  One  or  two  objects  in  our  final  sample  could  be  such 
contaminants,  but  because  we  have  three  different  realizations 
of  ten  objects  (3  x  10),  each  showing  similar  profiles  (details 
discussed  in  Section  4.2),  they  do  not  appear  to  affect  our  results. 

The  z  —  4  sample  has  415  objects,  while  the  z  —  5  sample 
has  265  objects.  In  Figure  2(a)  and  (b),  we  show  the  distribution 
of  the  FWHM  and  ellipticity,  e  =  (1  —  b/a),  measured  in 
each  of  the  two  samples  using  SExtractor  (Bertin  &  Arnouts 
1996).  We  further  constrained  our  samples  by  imposing  limits 
on  compactness  and  on  roundness  of  FWHM  ^  O'.' 3  and 
e  <  0.3.  Again,  this  is  to  minimize  the  probability  that  the 
z  —  4-5  candidates  are  significantly  dynamically  disturbed, 
and  to  maximize  the  probability  of  selecting  physically  similar 
objects.  Our  goal  is  to  find  the  visibly  most  symmetric,  least 
disturbed  systems  for  the  current  study.  This  sub-selection 
leaves  204  objects  in  the  z  —  4  sample  and  102  objects 
in  the  z  ~  5  sample.  Most  of  these  objects  are  faint,  and 
are  only  a  few  pixels  across  in  size  and,  hence,  have  larger 
uncertainties  in  their  measurements  of  FWHM  and  ellipticity. 
Therefore,  we  also  checked  our  objects  visually  to  eliminate 
any  possibility  of  our  selected  objects  being  contaminated  by 
unrelated  nearby  objects,  being  clearly  extended,  or  objects  with 
complex  morphologies. 

3.2.  z  —  6  Objects  (i'-Band  Dropouts) 

Yan  &  Windhorst  (2004b)  found  108  possible  5.5  ^  z  ^ 
6.5  candidates  in  the  HUDF  to  /«ab(^85o)  =  30.0  mag.  Bunker 
et  al.  (2004)  independently  found  the  brightest  54  of  these  108 
z  —  6  candidates  to  AB  =  28.5  mag.  Similarly,  deep  HST/ACS 
grism  spectra  of  the  HUDF  /'-band  dropouts  confirm  >  90%  of 
these  objects  at  AB  <  27.5  mag  to  be  at  z  ~  6  (Malhotra  et  al. 
2005,  Hathi  et  al.  2007).  Using  the  catalog  of  Yan  &  Windhorst 
(2004b),  we  extracted  108  postage  stamps,  each  51x51  pixels 
in  size,  from  the  HUDF  z'-band  image. 

As  for  the  z  —  4  and  z  ~  5  objects,  for  each  z  —  6  object  we 
measured  its  z'-band  FWHM  and  ellipticity  using  SExtractor. 
Figure  2(c)  shows  the  measured  ellipticity  versus  FWHM  for  all 
108  z  ~  6  candidates.  A  smaller  sample  of  67  objects  satisfies 
our  constraints  on  the  FWHM  and  ellipticity.  Further  visual 
inspection,  to  make  sure  that  our  sample  has  only  isolated, 
compact,  and  round  objects,  leaves  30  objects  in  our  z  —  6 
sample.  We  therefore  also  imposed  a  sample  size  of  30  objects 
on  the  two  lower  redshift  bins  after  visual  inspection. 


FWHM  (arcsec) 


Figure  2.  Ellipticity,  (1  —  b/a),  versus  object  FWHM.  for  all  z  —  4  (a),  z  —  5 
(b),  and  z  —  6  (c)  objects  selected  in  the  HUDF.  Measurements  were  performed 
in  i'-band  for  z  —  4  and  z  —  5  objects,  while  we  used  the  z'-band  for  z  —  6 
objects.  The  FWHM  of  a  stellar  image/PSF  is  ~3  pixels  or  Of! 09.  indicated 
by  the  leftmost  hatched  area  in  each  panel.  Objects  within  the  shaded  area 
meet  our  additional  selection  criteria  on  roundness  (e  z!  0.3)  and  compactness 
(FWHM  ^  0"3  or  10  pixels). 

The  results  in  this  paper  are  therefore  based  on  approximately 
(30/415)  ~  7%,  (30/265)  ~  11%,  and  (30/108)  ~  28%  of  the 
total  z  —  4,  5,  and  6  galaxy  populations. 

4.  RESUFTS 

4.1.  The  HUDF  Sky  Surface-Brightness  Level  and  its  rms 
Variation 

For  the  present  work,  it  is  critical  that  we  accurately  char¬ 
acterize  the  sky  background,  and  correctly  propagate  the  true 
ler  errors  due  to  the  subtraction  of  this  sky  background.  In  the 
following,  we  will  pursue  two  complementary  approaches  to 
determine  the  sky  surface  brightness,  and  compare  the  results. 
Here,  we  discuss  the  z'-band  measurements  in  detail. 

We  first  measured  the  sky  background  in  each  of  the 
415  z  —  4  object  stamps  (“local”  sky  measurements).  The 
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Figure  3.  Distribution  of  the  modal  sky-background  level  used  to  estimate 
the  la  uncertainty  in  that  level,  as  measured  in  the  415  z  —  4  object  stamps 
extracted  from  the  drizzled  HUDF  images  (a)  for  V-band,  (b)  for  /'-band,  and 
(c)  for  z'-band.  The  mean  (p)  and  the  sigma  (a)  of  the  best-fit  Gaussian  to  these 
distributions  are  also  shown  in  each  panel. 

Interactive  Data  Language  (IDL5)  procedure  SKY/MMM.pro6 
was  used  to  measure  the  sky  background.  This  procedure  is 
adapted  from  the  DADPHOT  (Stetson  1987)  routine  of  the  same 
name  and  works  as  follows.  First,  the  average  and  sigma  are 
obtained  from  the  sky  pixels.  Second,  these  values  are  used  to 
eliminate  outliers  with  a  low  probability.  Third,  the  values  are 
then  recomputed  and  the  process  is  repeated  up  to  20  iterations. 
If  there  is  a  contamination  due  to  an  object,  then  the  contami¬ 
nation  is  estimated  by  comparing  the  mean  and  median  of  the 
remaining  sky  pixels  to  get  the  true  sky  value.  The  output  of  this 
procedure  is  the  modal  sky  level  in  the  image. 

Figure  3(c)  shows  a  histogram  of  the  z'-band  modal  sky  values 
obtained  from  all  415  object  stamps  extracted  from  the  drizzled 
HUDF  images.  The  la  uncertainty  in  the  sky,  as^y,  determined 
from  a  Gaussian  fit  to  the  histogram,  is  2. 19  x  10-5  electrons  s-1 


5  IDL  website  http://www.ittvis.com/index.asp. 

6  Part  of  the  IDL  Astronomy  User’s  Library,  see: 
http://idlastro.gsfc.nasa.gov/homepage.html. 


in  the  z'-band.  The  sky-background  level  within  the  HUDF  was 
obtained  from  the  original  flat-fielded  ACS  images,  because 
the  final  co-added  HUDF  data  products  are  sky-subtracted.  The 
header  parameters  MDRIZSKY  and  EXPTIME  were  used  to 
obtain  the  actually  observed  sky  value.  MDRIZSKY  is  the 
sky  value  in  electrons  (e-)  computed  by  the  MultiDrizzle 
code  (Koekemoer  et  al.  2002),  while  EXPTIME  is  the  total 
exposure  time  for  the  image  in  seconds,  so  that  the  average  sky 
value  in  the  HUDF  has  the  units  of  e- s- 1 .  Figure  4(d)  shows  the 
histogram  of  the  sky  values  obtained  from  288  HUDF  z'-band 
flat-fielded  exposures.  The  average  value  of  the  sky  background, 
7sky,  is  0.0205  le“s_1.  That  sky  value  is  measured  from  the  flat- 
fielded  individual  ACS  images  with  pixel  sizes  of  0'.'05  pix-1 
and  hence,  in  the  following  calculations,  the  average  sky  value 
is  multiplied  by  a  factor  of  (0.030/0. 05)2  =  0.602  to  obtain  the 
corresponding  average  sky  value  for  the  HUDF  drizzled  pixel 
size  of  a.'030  pix-1 .  Using  these  values,  we  estimate  the  relative 
rms  random  sky-subtraction  error  as  follows: 


ss,ran 


°sky,ran 

-^sky 


2.19  x  10-5 
2.05  x  10-2  ■  0.602 


=  2.97  x  10-3. 


The  measured  average  sky-background  level  can  then  be  ex¬ 
pressed  as  the  z'-band  sky  surface  brightness  as  follows: 


/I,,  =  24.862  —  2.5  •  log 


/  0.0205-  0.602\ 
V  0.0302  ) 


=  22.577  ±  0.003  mag  arcsec  2, 


where  24.862  is  the  ACS/WFC  z'-band  AB  zero  point,  and 
0T030  pixel-1  is  the  drizzled  pixel  scale.  This  is  consistent 
with  the  values  obtained  by  extrapolating  the  on-orbit  BVI 
sky  surface  brightness  of  Windhorst  et  al.  (1994,  1998)  to  z', 
with  the  sky-background  estimates  from  the  ACS  Instrument 
Handbook  (Gonzaga  et  al.  2005),  and  with  the  colors  obtained 
by  convolving  the  filter  transmission  curves  with  the  solar 
spectrum.  Table  1  gives  the  measured  electron  detection  rate, 
surface  brightness,  and  colors  of  the  sky  background  with 
their  corresponding  errors  for  the  HUDF  BVi'z'  bands  as 
calculated  from  Figures  3  and  4.  The  contribution  of  the  zodiacal 
background  dominates  the  total  sky  background,  which  we  find 
to  be  only  ~10%  redder  in  (V-ir)  and  (i'-z!)  than  the  Sun. 
The  z'-band  surface  brightness  corresponding  to  the  1  a  sky- 
subtraction  uncertainty  is  therefore 

H,j  -  2.5  •  log(2ss,ran)  =  22.577  -  2.5  •  log(2.97  x  10-3) 

—  28.895  mag  arcsec-2. 

Next,  we  measure  the  sky  background  from  415  “blank” 
sky  stamps  (51  x  51  pixels)  distributed  throughout  the  HUDF 
(“global”  sky  measurements).  We  measure  the  sky  background 
using  the  same  IDL  algorithm  as  used  above. 

Figure  5(c)  shows  the  histogram  of  the  measured  z'-band 
modal  sky  values.  A  Gaussian  distribution  was  fit  to  this 
histogram,  giving  a  sky-sigma  of  2.00  x  10-5  e-  s-1.  The 
average  value  of  the  sky  remains  0.02051  e-  s-1  (Figure  4(d)). 
Using  these  values,  we  can  estimate  a  relative  rms  systematic 
sky-subtraction  error  as  follows: 


ss,sys  — 


^sky.sys 

/sky 


2.00  x  10-5 
2.05  x  10-2  •  0.602 


=  2.71  x  10-3. 


Since  the  z'-band  sky  surface  brightness  remains  22.577  mag 
arcsec-2,  this  gives  us  for  the  surface  brightness  corresponding 
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Figure  4.  The  actual  sky  values  measured  using  header  parameters  MDRIZSKY  and  EXPTIME  from  flat-fielded  HUDF  exposures:  (a)  for  the  R-band  using  112 
exposures,  (b)  for  the  V-band  using  112  exposures,  (c)  for  the  /'-band  using  288  exposures,  and  (d)  for  the  z'-band  using  288  exposures.  The  mean  (/x)  and  the  sigma 
(cr)  of  the  best-fit  Gaussian  to  these  distributions  are  shown  in  each  panel. 


Table  1 

Measured  Sky  Values  in  BVi'z'  (Filters)  for  the  HUDF 


HUDF 

filter 

Number  of 

exposures 

Mean  sky  value3 
(e-  s-1)  and  rms  errorb 

Sky  SBC 

(AB  mag  arcsec-2) 

Sky  color0 
(AB  mag) 

1<7  sky-subtraction 
error  (AB  mag  arcsec-2) 

B 

112 

0.015909  ±  0.000065 

23.664  ±0.003 

(B  -  V0sky  =  0.800 

29.85  ±  0.05 

V 

112 

0.070276  ±  0.000297 

22.864  ±  0.002 

(V  -  i')sky  =  0.222 

30.15  ±0.15 

i' 

288 

0.040075  ±  0.000088 

22.642  ±  0.002 

a’  ~  z'kky  =  0.065 

29.77  ±  0.20 

z' 

288 

0.020511  ±0.000047 

22.577  ±  0.003 

(V  -  z')sky  =  0.287 

28.95  ±  0.05 

Notes. 

3  From  Figure  4. 

b  Error  is  standard  deviation  of  the  mean  (crA/A). 

c  Sky  surface-brightness  values  and  colors  are  consistent  with  the  solar  colors  in  AB  mag  of  ( V  —  /')  =  0. 19,  ( V  —  z')  = 
0.21  and  (/'  —  z')  =  0.01  (except  for  bluest  color  ( B  —  V)),  and  is  dominated  by  the  zodiacal  background. 


to  the  la  sky-subtraction  uncertainty: 

/V  -  2.5  •  log(Ess,sys)  =  22.577  -  2.5  •  log(2.71  x  1(T3) 

=  28.995  mag  arcsec-2. 

From  these  two  complementary  approaches,  we  can  conclude 
that  all  surface  brightness  measurements  become  unreliable 
for  surface-brightness  levels  fainter  than  28.95  ±  0.05  mag 
arcsec-2  in  the  z'-band.  We  have  also  experimented  with  slightly 
larger  cutouts  (75  x  75  pixels  instead  of  51  x  51  pixels)  to  es¬ 
timate  the  sky-subtraction  error.  We  find  that  with  the  larger 
cutouts,  the  surface  brightness  corresponding  to  the  la  sky- 


subtraction  error  is  ~0. 1-0.2  mag  arcsec-2  fainter.  For  larger 
cutouts  we  expect  this  surface  brightness  to  be  ~0.4  mag 
fainter  but  we  find  about  0. 1-0.2  mag  fainter.  This  might  be 
because  of  residual  systematic  errors  in  the  HUDF  images. 
Therefore,  we  are  at  the  limit  of  accurately  measuring  this  sur¬ 
face  brightness  and,  hence,  we  will  here  quote  the  conservative 
brighter  limit  of  the  surface  brightness  corresponding  to  this 
la  sky-subtraction  error.  Expected  contributions  to  this  surface 
brightness  due  to  uncertainties  in  the  bias  level  determinations, 
which  correspond  to  ~0.001  counts  s-1  for  typical  HUDF  ex¬ 
posures  (A.  M.  Koekemoer,  private  communication),  are  less 
than  1%. 


COMPOSITE  SURFACE  BRIGHTNESS  PROFILES 


161 


No.  1,2008 


Sky  Value  (xlO-s  electrons/sec) 

Figure  5.  Distribution  of  the  modal  sky-background  level  used  to  estimate  the 
1  o  uncertainty  in  that  level,  as  measured  in  41 5  "blank”  5 1  x  51  pixel  sky  stamps 
extracted  from  the  drizzled  HUDF  images  (a)  for  V-band,  (b)  for  /'-band,  and 
(c)  for  r'-band.  The  mean  {// )  and  the  sigma  (cr)  of  the  best-fit  Gaussian  to  these 
distributions  are  shown  in  each  panel. 

Figure  5  clearly  shows  that  the  distribution  of  the  modal  sky 
values  is  not  as  symmetric  around  zero  as  in  Figure  3,  and  hence, 
the  use  of  a  “global”  sky  value  for  the  HUDF  is  not  as  reliable  as 
"local”  sky  measurements.  Therefore,  for  the  surface-brightness 
profiles  and  the  following  discussion,  we  will  adopt  the  local 
Ict  random  sky-subtraction  error  for  all  objects  in  our  study. 

The  average  modal  sky  values  and  their  1  a  errors  in  the  V- 
and  /'-bands  were  calculated  in  exactly  the  same  way  as  for 
the  z'-band,  as  shown  in  Figures  3-5.  The  resulting  BVi'z!  sky 
values  and  the  sky  surface-brightness  levels  are  all  given  in 
Table  1. 

4.2.  Composite  Images  and  Surface  Brightness  Profiles 

For  each  redshift  bin  (z  ~  4,  5,  and  6),  we  generated  three 
“stacked”  composite  images  from  subsets  of  ten  postage  stamps 
which  were  selected  as  follows.  After  placing  all  30  image 
stamps  per  redshift  bin  into  a  30  x  (51x51)  pixel  IDL  array, 
10  stamps  were  randomly  drawn  without  selecting  any  object 


more  than  once.  An  output  image  was  generated,  in  which 
the  values  at  each  pixel  are  the  average  of  the  corresponding 
pixels  in  the  ten  selected  input  stamps.  From  the  remaining  20 
stamps,  we  again  randomly  select  10,  from  which  we  generated 
a  second  composite  image,  after  which  the  final  10  images  were 
averaged  into  the  third  composite  image.  The  three  composite 
images  per  redshift  bin  are  therefore  independent  of  each  other. 
In  none  of  our  realizations  did  we  produce  composite  images  that 
were  essentially  unresolved.  Even  the  faintest  z  ~  4-6  galaxies 
are  clearly  resolved.  The  z  ~  4,  5,  and  6  objects  used  to  generate 
the  composite  images  have  an  apparent  magnitude  range  of 
approximately  27.5  ±1.0  AB  mag.  Because  the  magnitude 
range  is  relatively  small  and  the  S/N  per  pixel  is  low  even  in 
their  central  pixel,  we  have  given  all  objects  equal  weight.  To 
test  whether  this  range  in  magnitude  will  affect  our  stacks,  and 
hence  our  profiles,  we  created  three  stacks  depending  on  the 
apparent  magnitude,  i.e.,  one  stack  of  the  ten  brightest  objects 
in  the  sample,  a  second  stack  of  the  ten  next  brightest  objects 
in  the  sample  and  a  third  stack  of  the  ten  faintest  objects  in 
the  sample.  This  is  summarized  in  Figure  6(d).  We  found  that 
the  profiles  were  very  similar  except  that  the  profiles  of  the 
fainter  stacks  fall  off  more  quickly  at  larger  radius  compared 
to  the  profile  of  the  brightest  stack,  but  the  inner  profile  and 
the  deviation  in  the  profiles  are  clearly  visible  in  all  three 
stacks.  Therefore,  we  conclude  that  for  our  range  in  apparent 
magnitudes,  our  stacks/profiles  are  not  affected.  Perhaps  most 
surprisingly,  Figure  6(d)  shows  that  re  value  of  all  three  flux 
ranges  (~26. 0-27.0,  ~27. 0-28.0,  and  ~28. 0-29.0  mag)  are  all 
about  the  same  over  ~3-4  mag  in  flux,  so  the  primary  parameter 
that  distinguishes  the  brighter  from  the  fainter  z  ~  6  dropouts 
is  their  central  surface  brightness  (which  thus  also  varies  by 
~3-4  mag). 

We  used  the  IRAF7  procedure  ELLIPSE  to  fit  surface- 
brightness  profiles  shown  in  Figure  6  to  each  of  the  three 
independent  composite  images  per  redshift  bin.  We  also 
computed  a  mean  surface-brightness  profile  from  the  three  com¬ 
posite  surface-brightness  profiles  generated  from  the  three  in¬ 
dependent  composite  images  for  each  redshift  bin.  Figure  7 
shows  composite  images  for  z  ~  4,  5,  and  6  objects.  Here  each 
composite  image  is  a  stack  of  30  objects.  Figure  8  shows  the 
average  surface-brightness  profiles  for  each  of  the  redshift  in¬ 
tervals  z  ~  4,  5,  and  6.  The  thin  solid  curves  in  Figure  6  and  the 
dot-dash  curves  in  Figure  8  represent  the  observed  ACS  V-, 
and  z'-band  point  spread  functions  (PSFs),  while  the  horizontal 
dashed  lines  indicate  the  surface-brightness  level  corresponding 
to  the  Ict  sky-subtraction  error  in  each  of  the  HUDF  images  as 
discussed  in  Section  4.1.  It  is  important  to  note  that  we  scaled 
the  ACS  PSFs  to  match  the  surface  brightness  of  the  central  data 
point  in  our  mean  surface-brightness  profile,  to  determine  how 
extended  the  mean  surface-brightness  profile  is  with  respect  to 
the  PSFs. 

In  Figure  8,  we  fitted  all  possible  combinations  of  the  Sersic 
profiles  (convolved  with  the  ACS  PSF)  to  the  observed  profiles, 
and  using  x2  minimization  found  the  best  fits  for  galaxies  at 
z  —  4,  5,  and  6.  The  best-fit  Sersic  index  (n)  for  all  three 
profiles  (z  —  4,  5,  and  6)  is  n  <  2,  meaning  these  galaxies 
follow  mostly  exponential  disk-type  profiles  in  their  central 
regions.  We  find  that  the  observed  profiles  start  to  deviate 
from  the  best-fit  profiles  at  r  >  0'.'27,  somewhat  depending 


7  IRAF  (http://iraf.net)  is  distributed  by  the  National  Optical  Astronomy 
Observatory,  which  is  operated  by  the  Association  of  Universities  for  Research 
in  Astronomy,  Inc.,  under  cooperative  agreement  with  the  National  Science 
Foundation. 
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Semi-major  axis  [arcsec]  Semi-major  axis  [arcsec] 

Figure  6.  Composite  surface-brightness  profiles  for  three  independent  sets  of  ten  objects  at  (a)  z  —  4,  (b)  z  —  5,  and  (c)  z  —  6,  respectively.  The  thin  solid 
curve  represents  the  ACS  V-,  i'-,  and  "'-hand  PSFs,  respectively,  while  the  horizontal  dashed  line  indicates  the  surface-brightness  level  corresponding  to  the  1  cr 
sky-subtraction  error  in  the  HUDF  images.  The  vertical  dotted  line  marks  the  radius  at  which  the  profile  starts  to  deviate  significantly  from  the  extrapolation  of  the 
inner  r1'”  profile  observed  at  smaller  radii.  Note  that  at  z  —  6,  this  deviation  is  still  well  above  the  red  z'-band  PSF  halo  at  r  >  <y/30.  The  panel  (d)  shows  three 
z  —  6  composite  profiles  (each  with  a  set  of  ten  objects)  divided  by  apparent  magnitudes.  The  brightest  composite  profile  (dotted)  has  an  average  "'-band  magnitude 
of  ~26.8  mag.  The  next  brightest  composite  profile  (short  dash)  has  an  average  z'-band  magnitude  of  —27.9  mag,  and  the  faintest  composite  profile  (dot-dash)  has  an 
average  z'-band  magnitude  of  —28.9  mag. 


Figure  7.  Composite  images  for  z  ~  4  (left),  z  —  5  (center),  and  z  ~  6  (right)  objects.  Here  each  composite  image  is  a  stack  of  30  objects.  Each  stamp  is  1753  on  a 
side. 


on  the  redshift.  From  Figure  8,  we  also  see  that  in  each  of 
V  (z  —  4),  i'  ( z  —  5),  and  z!  ( z  —  6),  the  PSF  declines 
more  rapidly  with  radius  than  the  composite  radial  surface- 
brightness  profile  for  r  >  0'.'27.  It  is  therefore  unlikely  that  the 
observed  “breaks”  result  from  the  halos  and  structure  of  the 
ACS  PSFs.  Specifically,  at  z  —  6  the  most  significant  deviations 
in  the  light  profiles  are  seen  at  levels  1. 5-2.0  mag  above  the 
lcr  sky-subtraction  error,  and  well  above  the  PSF  wings.  Each 
of  the  mean  surface-brightness  profiles  displays  a  well-defined 
break,  the  radius  of  which  appears  to  change  somewhat  with 
redshift.  These  results  are  tabulated  in  Table  2.  The  vertical 


dotted  lines  (in  Figures  6  and  8)  mark  the  radius  at  which  the 
mean  surface-brightness  profiles  start  to  deviate  significantly 
from  the  extrapolation  of  the  r1/"  profile  observed  at  smaller 
radii. 

4.3.  Test  of  the  Stacking  Technique  on  Nearby  Galaxies 

To  test  the  general  validity  of  the  stacking  technique  itself 
on  a  local  galaxy  sample,  we  used  surface  photometry  from 
the  Nearby  Field  Galaxy  Survey  (NFGS;  Jansen  et  al.  2000a, 
2000b).  The  NFGS  sample  contains  196  nearby  galaxies, 
which  were  objectively  selected  from  the  CfA  redshift  catalog 
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Figure  8.  Mean  surface-brightness  profiles  with  a  best-fit  Sersic  profiles  for  30  composite  images  at  (a)  z  ~  4,  (b)  z  ~  5,  and  (c)  z  —  6,  respectively.  The  thin 
dot-dash  curve  represents  the  ACS  V-, and  ,7-band  PSFs,  respectively,  while  the  horizontal  dashed  line  indicates  the  surface-brightness  level  corresponding  to  the 
la  sky-subtraction  error  in  the  HUDF  images.  The  vertical  dotted  line  marks  the  radius  at  which  the  profile  starts  to  deviate  significantly  from  the  extrapolation  of  the 
inner  r1/"  profile  observed  at  smaller  radii.  The  n  and  rc  are  the  best-fit  Sersic  parameters. 

(CfA  I;  Davis  &  Peebles  (1983)  and  Huchra  et  al.  (1983) 
to  span  the  full  range  in  absolute  B  magnitude  present  in 
the  CfA  I  (—14.7  <  MB  <  -22.1  mag).  The  absolute 
magnitude  distribution  in  the  NFGS  sample  approximates  the 
local  galaxy  luminosity  function  (e.g.,  Marzke  et  al.  1994),  while 
the  distribution  over  Hubble  type  follows  the  changing  mix  of 
morphological  types  as  a  function  of  luminosity  in  the  local 
galaxy  population.  The  NFGS  sample  (as  detailed  in  Jansen 
et  al.  2000b)  minimizes  biases,  and  yields  a  sample  that,  with 
very  few  caveats,  is  representative  of  the  local  galaxy  population. 

As  part  of  the  NFGS,  U  BR  surface  photometry,  both  integrated 
(global)  and  nuclear  spectrophotometry,  as  well  as  internal 
kinematics  were  obtained  (see  Jansen  &  Kannappan  2001). 

Here,  we  will  concentrate  on  the  U -band  surface  photometry, 
since  it  is  closest  in  wavelength  to  the  rest-frame  wavelengths 
observed  at  z  —  4-6.  Although,  ideally,  we  would  want  a  filter 
further  into  the  UV,  Taylor  et  al.  (2007)  and  Windhorst  et  al. 

(2002)  show  that  for  the  majority  of  late-type  nearby  galaxies 
the  apparent  structure  of  galaxies  does  not  change  dramatically 
once  one  observes  shortward  of  the  B  aimer  break.  Early-type 
galaxies,  however,  are  a  clear  exception  to  this,  but  these  are 


Table  2 

Dynamical  Ages  for  z  —  4-6  Objects  in  the  HUDF 


Redshift 

z 

“Break”  radius3 
(arcsec) 

“Break”  radiusb 
(kpc) 

Dynamical  agec 

(Tdyn) 

4 

0.35 

2.5 

0.09-0.29  Gyr 

5 

0.31 

2.0 

0.07-0.21  Gyr 

6 

0.27 

1.6 

0.05-0.15  Gyr 

Notes. 

“From  composite  surface-brightness  profiles  (Figures  6  and  8). 
bRadius  in  kpc  corresponding  to  radius  in  arcsec  at  given  redshift. 
cIf  “break  radius”  interpreted  as  an  indicator  of  dynamical  age. 

not  believed  to  dominate  the  galaxy  population  at  z  ~  4-6,  as 
discussed  before. 

Figure  9  shows  stacked  profiles  for  relatively  luminous  early-, 
spiral-,  and  late-type  galaxies  drawn  from  the  NFGS.  Vertical 
dotted  lines  indicate  the  half-light  radii  and  their  intersection 
with  the  profiles,  the  surface  brightness  at  that  radius.  Dashed 
lines  indicate  exponential  fits  to  the  outer  portion  of  each 
profile.  Figure  9  also  shows  that  co-adding  profiles  for  disparate 
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Figure  9.  Stacked  radial  surface-brightness  profiles  for  relatively  luminous  early-,  spiral-,  and  late-type  nearby  galaxies  drawn  from  the  Nearby  Field  Galaxy  Survey 
(Jansen  et  al.  2000a,  2000b).  The  vertical  dotted  line  indicates  the  half-light  radius,  while  the  dashed  line  represents  an  exponential  fit  to  the  outer  portions  of  each 
composite  profile.  Co-adding  profiles  for  disparate  morphological  types  and  for  spiral  galaxies  with  a  range  in  bulge-to-disk  ratios  can  produce  breaks  in  the  composite 
profile.  No  significant  breaks  are  seen  in  the  outer  light  profiles,  when  the  profiles  of  either  early-type  galaxies  (E,  SO)  or  late-type  galaxies  (Sd-Irr)  are  co-added. 


morphological  types  and  for  mid-type  spiral  galaxies  with 
a  range  in  bulge-to-disk  ratios  can  produce  breaks  in  the 
composite  profile.  No  such  breaks  are  seen  when  the  profiles  of 
either  early-type  galaxies  (E,  SO)  or  late-type  galaxies  (Sd-Irr) 
are  co-added.  This  figure  shows  that  if  galaxies  at  z  —  4-6  had 
similar  morphological  types  as  local  galaxies,  then  it  would  be 
possible  to  produce  a  break  in  the  profiles  (as  shown  in  Figures  6 
and  8),  merely  by  mixing  different  types  of  galaxies.  We  do  not 
believe  that  the  galaxy  populations  at  z  —  4-6  morphologically 
resemble  those  at  low  redshift.  Hence,  for  primarily  late-type 
galaxies,  which  dominate  the  faint  blue  galaxy  population  at 
AB  >  24  mag  (Driver  et  al.  1 998),  and  which  likely  dominate  the 
fainter  end  of  the  luminosity  function  at  z  —  4-6  that  we  sample 
here  (Yan  &  Windhorst  2004a,  2004b),  the  image  stacking  is 
likely  a  valid  exercise. 

The  primary  goal  of  this  section  was  to  show  that  the  profile 
stacking  technique  is  valid  and  can  be  used  to  get  meaningful 
surface-brightness  profiles.  We  are  not  comparing  our  nearby 
sample  with  galaxies  at  z  —  4-6.  These  nearby  galaxies  are 
unlikely  to  be  local  analogs  of  high-redshift  galaxies.  If  we 
apply  surface  brightness  dimming  to  UV  light  profiles  of  these 
nearby  galaxies,  they  would  be  mostly  invisible  to  the  HST, 
and  in  some  cases  visible  to  the  James  Webb  Space  Telescope 
(. JWST)  in  long  integration  (see,  e.g.,  Windhorst  et  al.  2006). 


This  is  another  way  of  saying  that  the  z  —  4-6  objects  are  truly 
different  from  the  z  —  0  objects. 

5.  DISCUSSIONS 

Figure  8  shows  that  the  mean  surface-brightness  pro¬ 
files  deviate  significantly  from  an  inner  r  profile  at  radii 
r  >  0'.'27-0"35,  depending  somewhat  on  the  redshift  bin.  These 
deviations  appear  real,  with  the  break/point  of  departure  located 
>1.5-2  mag  above  the  la  sky-subtraction  error  and  above  the 
PSF  wings.  In  the  following,  we  discuss  several  possible  ex¬ 
planations  for  the  observed  shapes  of  our  composite  surface- 
brightness  profiles. 

5.1.  Galaxies  with  Different  Morphologies 

Our  test  on  nearby  galaxies  (Figure  9)  shows  that  if  we 
stack  many  galaxies  with  different  morphologies  (early-type, 
late-type,  or  spiral  galaxies),  it  is  possible  to  get  a  slope 
change  (“break”)  in  the  average  surface-brightness  profile. 
Ravindranath  et  al.  (2006)  find  that  40%  of  the  brighter 
FBGs  at  2.5  <  z  <  5  have  light  profiles  close  to  expo¬ 
nential,  as  seen  for  disk  galaxies,  and  only  ~30%  have  high 
n,  as  seen  in  nearby  spheroids.  They  also  find  a  significant 
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fraction  (~30%)  of  galaxies  with  light  profiles  shallower  than 
exponential,  which  appear  to  have  multiple  cores  or  disturbed 
morphologies,  suggestive  of  close  pairs  or  on-going  galaxy 
mergers.  Therefore,  if  galaxies  at  z  —  4-6  have  a  variety  of 
morphological  types,  then  the  shape  of  the  average  surface- 
brightness  profile  that  we  see  may  be  due  to  the  stacking  of 
different  types  of  galaxies.  Therefore,  we  find  that  the  exponen¬ 
tial  and  the  flatter  profiles  found  by  Ravindranath  et  al.  (2006) 
for  galaxies  at  2.5  <  z  <  5  also  apply  to  higher  redshifts 
(z  >  5). 

Also,  we  believe  that  it  is  more  likely  that  the  high  redshift, 
faint  galaxy  population  consists  primarily  of  small  galaxies  with 
late-type  morphologies  and  with  sub-L*  luminosities,  as  seen  at 
z  —  2-3  (Driver  et  al.  1995, 1998).  So  if  the  z  —  4-6  population 
consists  of  such  a  late-type  galaxy  population,  then  the  slope 
change  in  the  light  profiles  is  likely  not  the  result  of  co-adding 
images  of  objects  with  disparate  morphological  types. 

5.2.  Central  Star  Formation/Starburst 

HST  optical  images  of  galaxies  at  z  —  4-6  sample  their 
rest-frame  UV  (~1200  A),  where  the  contribution  from  the 
actively  star-forming  regions  (very  young,  massive  stars)  dom¬ 
inates  the  UV  light.  Hathi  et  al.  (2007)  have  shown  that  galax¬ 
ies  at  z  —  5-6  are  high-redshift  starbursts  and  these  galax¬ 
ies  have  similar  starburst  intensity  limit  as  local  starbursting 
galaxies.  Therefore,  it  is  possible  that  galaxies  at  z  —  4-6 
have  centrally  concentrated  star  formation  or  starburst.  This 
possibility  is  based  on  three  key  assumptions:  (1)  most  of  the 
galaxies  at  z  —  4,  5,  and  6  are  intrinsically  late-type  galaxies 
(Driver  et  al.  1998;  Steidel  et  al.  1999);  (2)  the  spectral  en¬ 
ergy  distribution  (SED)  of  these  galaxies  at  z  —  4,  5,  and  6  is 
dominated  by  early-A-  to  late-O-type  stars,  respectively;  and 
(3)  there  are  no  old  stars  with  ages  at  z  —  4-6  greater  than 
2-1  Gyr  in  WMAP  cosmology,  respectively. 

Hunter  &  Elmegreen  (2006)  studied  azimuthally  averaged 
surface  photometry  profiles  for  large  sample  of  nearby  irregular 
galaxies.  They  found  that  some  galaxies  have  double  exponen¬ 
tials  that  are  steeper  (and  bluer)  in  the  inner  parts  compared  to 
outer  parts  of  the  galaxy.  They  discussed  that  this  type  of  be¬ 
havior  is  expected  in  galaxies  where  the  centrally  concentrated 
star  formation  or  starburst  steepens  the  surface-brightness  pro¬ 
files  in  the  center.  If  that  is  the  case,  then  one  might  expect  a 
better  correlation  between  the  break  in  the  surface-brightness 
profiles  and  changes  in  color  profiles.  Unfortunately,  for  our 
sample  of  galaxies  at  z  —  4-6,  we  do  not  have  high-resolution 
rest-frame  UBV  color  information.  The  objects  are  generally 
too  faint  for  the  Spitzer  Space  Telescope,  and  hence  we  cannot 
confirm  or  reject  this  possibility  for  the  shape  of  our  composite 
surface-brightness  profiles. 

5.3.  Limits  to  Dynamical  Ages  for  z  —  4,  5,  and  6  Objects 

The  average  compact  z  —  4-6  galaxy  is  clearly  extended 
with  respect  to  the  ACS  PSFs  (Figure  8),  and  is  best  fit  by  an 
exponential  profile  ( n  <  2)  out  to  a  radius  of  about  r  ~  0'.'35, 
0731,  and  0727  at  z  —  4,  5,  and  6,  respectively.  The  apparent 
progression  with  redshift  is  noteworthy.  The  radius  at  which  the 
profile  starts  to  deviate  from  r1!"  (in  this  case  at  radius  r  >  0737- 
0727)  may  put  an  important  constraint  on  the  dynamical  time 
scale  of  the  system,  as  discussed  in  Section  1 .  If  this  argument 
is  valid,  then  we  can  estimate  limits  to  the  dynamical  ages  of 
z  —  4,  5,  and  6  galaxies  as  follows. 

In  WMAP  cosmology,  a  radius  of  r  >  0735  at  z  —  4 
corresponds  to  r  >  2.5  kpc.  The  dynamical  time  scale  (e.g.. 


Binney  &  Tremaine  1987),  r^y,,,  goes  as  Tdyn  =  Cr3/2/V G  M, 
where  the  constant  C  =  n / 2.  For  a  typical  dwarf  galaxy  mass 
range  of  ~109-108Mo  inside  r  =  2.5  kpc,  we  infer  that  the 
limits  to  the  dynamical  age  would  be  rdyn  —  90-290  Myr,  which 
is  the  lifespan  expected  for  a  late-type  B-star.  This  means  that 
the  last  major  merger  that  affected  this  surface-brightness  profile 
and  that  triggered  its  associated  starburst  may  have  occurred 
~0.20  Gyr  before  z  ~  4 — assuming  that  the  star  formation 
was  not  spontaneous,  but  associated  with  some  accretion  or  a 
merging  event. 

Table  2  shows  the  break  radius  and  inferred  limits  to  dynam¬ 
ical  ages  for  the  z  ~  4-6  objects.  At  z  —  5,  we  find  that  the 
limits  to  dynamical  age  at  the  break  radius  would  be  rdyn  —  70- 
210  Myr,  which  is  the  lifespan  expected  for  a  mid-B-star,  while 
at  z  —  6,  tdyn  —  50-150  Myr,  which  is  the  lifespan  expected 
for  a  late-O-early-B-star.  This  means  that  the  last  major  merger 
that  affected  these  surface-brightness  profiles  at  z  —  5  and  6  and 
that  triggered  its  associated  starburst  may  have  occurred  ~0. 14 
and  ~0.10  Gyr  before  z  —  5  and  6,  respectively. 

The  dynamical  time  is  a  lower  limit  to  the  actual  time 
available,  since  it  assumes  matter  starts  from  rest.  Any  angular 
momentum  at  the  start  will  increase  the  available  time.  The  best- 
fit  SED  age  from  the  GOODS  HST  and  Spitzer  photometry  on 
some  of  the  brighter  of  these  objects — using  Bruzual  &  Chariot 
(2003)  templates — is  in  the  range  of  about  ~  150-650  Myr  (Yan 
et  al.  2005;  Eyles  et  al.  2005,  2007),  the  lower  end  of  which 
is  consistent  with  our  limits  to  their  dynamical  age  estimates, 
while  the  somewhat  larger  SED  ages  could  also  be  affected  by 
the  onset  of  the  AGB  in  the  stellar  population  increasing  the 
observed  Spitzer  fluxes  and  hence  possibly  overestimating  ages 
(Maraston  2005).  Our  age  estimates  for  z  —  4-6  are  consistent 
with  the  trend  of  SED  ages  suggested  for  z  —  7  (Labbe 
et  al.  2006).  It  is  noteworthy  that,  given  the  uncertainties,  the 
two  independent  age  estimates  are  consistent.  If  our  limits  to 
dynamical  age  estimates  for  the  image  stacks  are  thus  valid, 
they  are  consistent  with  the  SED  ages,  and  point  to  a  consistent 
young  age  for  these  objects. 

Furthermore,  the  presence  of  young,  massive  late-O-early- 
B-stars  at  z  —  6  has  implications  for  the  reionization  of  the 
Universe.  From  observations  of  the  appearance  of  complete 
Gunn-Peterson  troughs  in  the  spectra  of  z  >  5.8  quasars 
(Fan  et  al.  2006),  we  know  that  the  epoch  of  reionization 
had  ended  by  z  —  6.  From  the  steep  (a  =  —1.8)  faint-end 
slope  of  the  luminosity  function  of  z  ~  6  galaxies,  Yan  & 
Windhorst  (2004a,  2004b)  concluded  that  dwarf  galaxies,  and 
not  quasars,  likely  finished  reionization  by  z  ~  6.  Should  the 
present  interpretation  of  their  light  profiles  be  correct,  then  it 
would  appear  to  add  support  to  this  picture,  in  the  sense  that 
such  objects  are  dominated  by  B-stars  and  did  not  start  their 
most  recent  major  starburst  long  before  z  —  6. 

6.  SUMMARY 

We  used  the  stacked  HUDF  images  to  analyze  the  average 
surface-brightness  profiles  of  z  ~  4-6  galaxies.  Our  analysis 
shows  that  even  the  faintest  galaxies  at  z  —  4-6  are  resolved. 
This  may  have  implications  on  the  stellar  density  and  its  relation 
to  the  stellar  density  in  present-day  galaxies.  We  also  find  that 
the  average  surface-brightness  profiles  display  breaks  at  a  radius 
that  progresses  toward  lower  redshift  from  r  ~  0727  (1.6  kpc) 
at  z  —  6  to  r  ~  0735  (2.5  kpc)  at  z  ~  4. 

The  shape  of  the  radial  surface-brightness  profile  that  we  ob¬ 
serve  could  result  from  a  mixture  of  different  morphological 


166 


HATHI  ET  AL. 


Vol.  135 


types  of  galaxies,  if  they  exist  at  z  —  4-6,  because  we  can  pro¬ 
duce  similar  breaks  in  the  surface-brightness  profiles  when  we 
mix  different  types  of  nearby  galaxies.  Alternatively,  if  these 
galaxies  are  dominated  by  a  central  starburst  then  they  could 
show  such  double  exponential-type  profiles,  as  discussed  by 
Hunter  &  Elmegreen  (2006).  In  a  third  scenario,  if  the  galax¬ 
ies  at  z  —  4-6  are  truly  young  and  mostly  late  type,  the  outer 
profiles  seen  in  our  mean  radial  surface-brightness  profiles  at 
z  —  4-6  bear  the  imprint  of  the  hierarchical  build-up  pro¬ 
cess,  and  are  still  dominated  by  infalling  material,  which  is 
not  detectable  in  the  individual  HUDF  images  of  these  faint 
objects.  We  have  estimated  limits  to  dynamical  ages  from  the 
break  radius  at  z  —  4,  5,  and  6,  very  roughly  as  ~0.20,  0.14, 
and  0.10  Gyr,  respectively,  and  those  ages  are  similar  to  the 
SED  ages  inferred  at  z  —  4-6  (Yan  et  al.  2005;  Eyles  et  al. 
2005,  2007),  and  consistent  with  SED  ages  suggested  for  z  —  7 
(Labbe  et  al.  2006).  Hence,  at  z  —  4,  5,  and  6,  the  last  ma¬ 
jor  merger  that  affected  the  surface  brightness  profiles  that  we 
observe,  and  that  triggered  the  observed  star-burst,  may  have 
occurred  respectively  ~0.20,  0.14,  and  0.10  Gyr  earlier,  or  very 
approximately  at  z  —  4.5,  5.5,  and  6.5.  This  would  be  consistent 
with  the  hierarchical  assembly  of  galaxies  and  with  the  end  of 
reionization,  since  it  would  imply  that  from  z  —  Atoz  —  6  the 
SEDs  become  progressively  more  dominated  by  late-B-late-0 
stars.  This  implies  that  the  sub-L*  (i.e.  dwarf)  galaxies  may 
have  produced  sufficient  numbers  of  energetic  UV  photons  to 
complete  the  reionization  process  by  z  —  6,  as  Yan  &  Windhorst 
(2004a,  2004b)  suggested.  It  will  be  imperative  to  study  with 
future  instruments  like  //Y77WFC3  and  JWST  (Windhorst  et 
al.  2006,  2007)  whether  the  dominant  stellar  population  indeed 
changes  from  late-O-early-B  at  z  ~  6  (i.e.  capable  of  reioniz¬ 
ing)  to  mid-  to  late-B  at  z  —  4-5  (i.e.  capable  of  maintaining 
reionization),  and  to  what  extent  the  intrinsic  sizes  of  these  faint 
objects  will  ultimately  limit  deep  JWST  surveys. 

This  work  was  partially  supported  by  HST  grants  AR  10298 
and  GO  9780  from  the  Space  Telescope  Science  Institute,  which 
is  operated  by  AURA  under  NASA  contract  NAS  5-26555.  The 
authors  thank  Deidre  Hunter,  Alan  Dressier,  Henry  Ferguson, 
Anton  Koekemoer,  and  Robert  Morgan  for  their  helpful  discus¬ 
sions.  R.A.W.  acknowledges  support  from  NASA  JWST  Inter¬ 
disciplinary  Scientist  grant  NAG5- 12460  from  GSFC,  which 
supported  an  investigation  of  the  implications  of  this  work  for 
JWST.  We  specially  thank  our  referee,  Dr.  Patrick  McCarthy, 
for  his  helpful  comments  that  have  improved  this  paper. 
Facilities:  HST(ACS) 

REFERENCES 

Beckwith,  S„  et  al.  2006,  AJ,  132,  1729 
Bertin,  E.,  &  Arnouts,  S.  1996,  A&AS,  1 17,  393 

Binney,  J.  J.,  &  Tremaine,  S.  1987,  Galactic  Dynamics  (Princeton,  NJ:  Princeton 
Univ.  Press) 

Bolzonella,  M„  Miralles,  J.  M„  &  Pello,  R.  2000,  A&A,  363,  476 
Bouwens,  R.,  et  al.  2004,  ApJ,  606,  L25 

Bouwens,  R.  J.,  Illingworth,  G.  D.,  Blakeslee,  J.  P.,  &  Franx,  M.  2006,  ApJ, 
653,  53 

Bouwens,  R.  J.,  Illingworth,  G.  D..  Franx,  M.,  &  Ford,  H.  2007,  ApJ  in  press 
(arXiv:  0707.2080) 

Brandt,  W.  N.,  Hornschemeier,  A.  E.,  Schneider,  D.  P.,  Alexander,  D.  M.,  Bauer, 
F.  E„  Garmire,  G.  P.,  &  Vignali,  C.  2001,  ApJ,  558,  L5 
Bruzual,  G.,  &  Chariot,  S.  2003,  MNRAS,  344,  1000 

Bunker,  A.  J.,  Stanway,  E.  R.,  Ellis,  R.  S.,  &  McMahon,  R.  G.  2004,  MNRAS, 
355, 374 

Cowie,  L.  L„  Songaila,  A.,  Hu,  E.  M„  &  Cohen,  J.  G.  1996,  AJ,  112,  839 
Davis,  M„  &  Peebles,  P.  J.  M.  1983,  ApJ,  267,  465 


Dow-Hygelund,  C.,  et  al.  2007,  ApJ,  660,  47 

Driver,  S.  R,  Fernandez-Soto,  A.,  Couch,  W.  J.,  Odewahn,  S.  C.,  Windhorst, 

R.  A.,  Phillips,  S.,  Lanzetta,  K„  &  Yahil,  A.  1998,  ApJ,  496,  L93 
Driver,  S.  P„  Windhorst,  R.  A.,  &  Griffiths,  R.  E.  1995,  ApJ,  453,  48 

Eyles,  L.  P.,  Bunker,  A.  J.,  Ellis,  R.  S.,  Lacy,  M.,  Stanway,  E.  R.,  Stark,  D.,  & 
Chiu,  K.  2007,  MNRAS,  374,  910 

Eyles,  L.  P.,  Bunker,  A.  J.,  Stanway,  E.  R.,  Lacy,  M„  Ellis,  R.  S.,  &  Doherty,  M. 

2005,  MNRAS,  364,  443 
Fan,  X.,  et  al.  2006,  AJ,  132,  1 17 

Georgakakis,  A.,  Hopkins,  A.  M.,  Sullivan,  M.,  Afonso,  J.,  Georgantopoulos, 
I.  JMobasher,  B„  &  Cram,  L.  E.  2003,  MNRAS,  345,  939 
Giavalisco,  M.,  et  al.  2004,  ApJ,  600,  L103 

Gonzaga,  S.,  et  al.  2005,  in  ACS  Instrument  Handbook,  Version  6.0,  (Baltimore, 
MD:  STScI) 

Guzman,  R.,  Gallego,  J.,  Koo,  D.  C.,  Phillips,  A.  C.,  Lowenthal,  J.  D.,  Faber, 

S.  M„  Illingworth,  G.  D„  &  Vogt,  N.  P.  1997,  ApJ,  489,  559 

Hathi,  N.  P.,  Malhotra,  S.,  &  Rhoads,  J.  E.  2007,  ApJ,  in  press 
(arXiv:  0709.0520) 

Heavens,  A.,  Panter,  B.,  Jimenez,  R.,  &  Dunlop,  J.  2004,  Nature,  428,  625 
Hu,  E.  M.,  Cowie,  L.  L.,  McMahon,  R.  G.,  Capak,  P.,  Iwamuro,  F.,  Kneib,  J.-P.. 

Maihara,  T.,  &  Motohara,  K.  2002,  ApJ,  568,  L75 
Huchra,  J.  P.,  Davis,  M.,  Latham,  D..  &  Tonry,  J.  1983,  ApIS,  52,  89 
Hunter,  D.  A.,  &  Elmegreen,  B.  G.  2006,  ApJS,  162,  49 
Jansen,  R.  A.,  Fabricant,  D„  Franx,  M„  &  Caldwell,  N.  2000a,  ApJS, 
126,  331 

Jansen,  R.  A.,  Franx,  M.,  Fabricant,  D.,  &  Caldwell,  N.  2000b,  ApJS, 
126,  271 

Jansen,  R.  A.,  &  Kannappan,  S.  J.  2001,  Ap&SS,  276,  1151 
Kodaira,  K„  et  al.  2003,  PASJ,  55,  L17 
Kodama,  T„  et  al.  2004,  MNRAS,  350,  1005 

Koekemoer,  A.  M.,  Fruchter,  A.  S.,  Hook,  R.  N.,  &  Hack,  W.  2002,  ed. 
S.  Arribas,  A.  Koekemoer,  &  B.  Whitmore  The  2002  HST  Calibration 
Workshop  (Baltimore,  MD:  STScI),  350,  337 
Kormendy,  J.  1977,  ApJ,  218,  333 

Kurk,  J.  D.,  Cimatti,  A.,  di  Serego,  A.  S.,  Vernet,  J.,  Daddi,  E.,  Ferrara,  A.,  & 
Ciardi,  B.  2004,  A&A,  422,  L13 

Labbe,  I.,  Bouwens,  R.,  Illingworth,  G.  D.,  &  Franx,  M.  2006,  ApJ,  649,  L67 

Lynden-Bell,  D.  1967,  MNRAS,  136,  101 

Malhotra,  S.,  et  al.  2005,  ApJ,  626,  666 

Maraston,  C.  2005,  MNRAS,  362,  799 

Marzke,  R.  O.,  Huchra,  J.  P„  &  Geller,  M.  J.  1994,  ApJ,  428,  43 

McCarthy,  P.  J.  2004,  BAAS,  36,  1555 

Nandra,  K.,  Mushotzky,  R.  F.,  Arnaud,  K.,  Steidel,  C.  C.,  Adelberger,  K.  L., 
Gardner,  J.  P,  Teplitz,  H.  I.,  &  Windhorst,  R.  A.  2002,  ApJ,  576,  625 
Oke,  J.  B„  &  Gunn,  J.  E.  1983,  ApJ,  266,  713 

Panter,  B.,  Jimenez,  R.,  Heavens,  A.  F.,  &  Chariot,  S.  2007,  MNRAS, 
378, 1550 

Pascarelle,  S.  M.,  Windhorst,  R.  A.,  Keel,  W.  C.,  &  Odewahn,  S.  C.  1996,  Nature, 
383,  45 

Ravindranath,  S.,  et  al.  2006,  ApJ,  652,  963 
Rhoads,  J.  E„  et  al.  2004,  ApJ,  61 1,  59 
Spergel,  D.  N„  et  al.  2007,  ApJS,  170,  377 

Steidel,  C.  C.,  Adelberger,  K.  L.,  Giavalisco,  M.,  Dickinson,  M.,  &  Pettini,  M. 
1999,  ApJ,  519,  1 

Stern,  D.,  Yost,  S.  A.,  Eckart,  M.  E.,  Harrison,  F.  A.,  Helfand,  D.  J.,  Djorgovski, 
S.  G.,  Malhotra,  S.,  &  Rhoads,  J.  E.  2005,  ApJ,  619,  12 
Stetson,  P.  B.  1987,  PASP,  99,  191 
Taniguchi,  Y„  et  al.  2005,  PASJ,  57,  165 

Taylor,  V.  A.,  Conselice,  C.  J.,  Windhorst,  R.  A.,  &  Jansen,  R.  A.  2007,  ApJ, 
659,  162 

van  Albada,  T.  S.  1982,  MNRAS,  201,  939 
van  Dokkum,  P.  G.,  et  al.  2003,  ApJ,  587,  L83 
van  Dokkum,  P.  G.,  et  al.  2004,  ApJ,  611,  703 
Vanzella,  E„  et  al.  2006,  A&A,  454,  423 

White,  R.  L.,  Helfand,  D.  J.,  Becker,  R.  H.,  Glikman,  E.,  &  de  Vries,  W. 
2007,  ApJ  ,  654,  99 

Williams,  R.  E„  et  al.  1996,  AJ,  1 12,  1335 

Windhorst,  R.  A.,  Cohen,  S.  H.,  Jansen,  R.  A.,  Conselice,  C.,  &  Yan,  H. 

2006,  NewAR,  50, 1 13 

Windhorst,  R.  A.,  Gordon,  J.  M.,  Pascarelle,  S.  M.,  Schmidtke.  P  C.,  Keel,  W. 

C.,  Burkey,  J.  M„  &  Dunlop,  J.  S.  1994,  ApJ,  435,  577 
Windhorst,  R.  A.,  Hathi,  N.  P.,  Cohen,  S.  H.,  &  Jansen,  R.  A.  2007,  Advances 
in  Space  Research,  in  press  (arXiv  astro-ph/0703171) 

Windhorst,  R.  A.,  Keel,  W.  C.,  &  Pascarelle,  S.  M.  1998,  ApJ,  494,  27 

Windhorst,  R.  A.,  et  al.  2002,  ApJS,  143,  113 

Yan,  H„  et  al.  2005,  ApJ,  634,  109 

Yan,  H„  &  Windhorst,  R.  2004a,  ApJ,  600,  LI 

Yan,  H„  &  Windhorst,  R.  2004b,  ApJ,  612,  L93 

York,  D.  G.,  et  al.  2000,  AJ,  120,  1579 

Zibetti,  S.,  White,  S.  D.  M„  &  Brinkmann,  J.  2004,  MNRAS,  347,  556 


